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After  two  previous  reports  devoted  to  bibliographic  study  (interim  report  #  1)  and 
thermodynamic  calculations  (interim  report  #  2),  this  final  report  examines  the 
thermochemistry  and  reactivity  of  gaseous  mixtures  formed  above  ionic  liquids  (ILs)  when 
mixed  with  white  fuming  nitric  acid  (WFNA).  After  a  general  introduction  on  the  interesting 
results  recently  observed  in  this  field  by  Chambreau  and  Vaghjiani  et  al.  (2008)  at  the  Air 
Force  Research  Laboratory  (AFRL),  Edwards  AFB,  this  report  will  discuss  the 
thermochemistry  of  imidazoles,  triazoles,  tetrazoles,  and  nitro  compounds,  which  are 
energetic  compounds  for  which  relatively  little  thermodynamic  data  is  available.  Then,  the 
reactivity  of  the  gas  phase  above  the  reacting  ILs-WFNA  liquid  mixtures  will  be  examined  by 
performing  calculations  at  thermodynamic  equilibrium  as  described  in  interim  report  #  2. 
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Lastly,  a  detailed  chemical  kinetics  model  will  be  proposed  to  tentatively  explain  the  gas 
phase  processes  that  lead  to  ignition  in  such  systems.  To  our  knowledge,  this  model 
represents  the  first  ever  attempt  to  predict  hypergolic  ignition  delay  times  in  ILs-WFNA 
mixtures.  Our  preliminary  sensitivity  analysis  study  of  the  model  provides  information  on 
reaction  steps  that  may  be  important  in  determining  the  ignition  delay.  Future  efforts  to  refine 
the  present  reaction  mechanism  data  set,  together  with  incorporation  of  new  experimental 
kinetics  information  as  it  becomes  available,  as  well  as  inclusion  of  fluid  dynamical  effects 
will  significantly  enhance  the  fidelity  of  the  model  to  quantitatively  predict  observed  ignition 
delays  in  ILs-WFNA  hypergol  mixtures  under  typical  laboratory  conditions. 


1 .  Introduction 

The  beginning  of  the  21st  century  is  particularly  interesting  from  an  energetics  point  of  view. 
On  one  hand,  the  increase  in  pollution  and  greenhouse  effects  can  be  correlated  to  the  rapid 
increase  in  combustion  of  fossil  fuels.  On  the  other  hand,  serious  depletions  of  known 
reservoirs  and  scarcity  of  easily  extractable  new  fossil  fuel  sources  has  led  to  the  search  for 
new  high-energy  density  and  environmentally  benign  fuels.  These  new  fuels  are  not  totally 
“green”  but  they  address,  at  least  partially,  some  of  the  drawbacks  of  current  fuels  and 
combustibles.  It  is  the  case  for  gasoline  and  diesel  that  bio-ethanol  and  methyl  esters  of 
vegetable  oils,  respectively,  will  be  the  likely  replacement  candidates.  This  ecological 
concern  is  not  limited  to  automobile  fuels  but  extends  to  all  fuels  and  combustibles,  including 
space  fuels,  although,  the  quantities  needed  for  space  fuels  are  somewhat  marginal  compared 
to  ground  transportation  fuels.  One  of  the  most  utilized  bipropellant  combinations  for 
propulsion  of  satellites  and  for  interplanetary  missions  that  fall  under  this  scope  is 
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monomethylhydrazine/nitrogen  tetroxide  (MMH/NTO).  These  propellants  have  the 
advantage  of  being  earth-storable  liquids  (non-cry ogenic);  however,  their  vapors  are  highly 
toxic  and  carcinogenic,  and  their  time-weighted  average  (TWA)  exposure  and  threshold  limit 
value  (TLV)  are  rather  low  for  these  compounds.  Moreover,  because  MMH  has  a  positive 
enthalpy  of  formation,  it  is  able  to  sustain  a  decomposition  flame,  and  therefore  can  be 
hazardous  even  in  the  absence  of  an  oxidizer.  Furthermore,  MMH  compatibility  is  poor  with 
many  materials  on  which  it  can  rapidly  decompose  exothermically  upon  contact.  Similarly, 
NTO  is  corrosive  to  many  materials.  It  is  clear  from  the  above  facts  that  there  are  numerous 
valid  arguments  to  search  for  suitable  replacements  for  MMH.  Finding  a  replacement  for 
NTO  is  also  of  interest  but  probably  less  crucial.  For  instance,  NTO  could  be  replaced  by 
nitric  acid  (or  by  any  other  suitable  oxidant)  when  necessary.  Among  some  of  the  alternative 
replacements  proposed  for  MMH  are  ionic  liquids  which  exhibit  much  lower  vapor  pressures 
than  MMH,  and  have  better  TWA  and  TLV  characteristics. 

However,  the  main  interesting  feature  of  the  MMH/NTO  bipropellant  combination  is  that 
MMH  and  NTO  form  mixtures  that  autoignite  at  room  temperature  and  at  low  pressures. 
Understanding  this  autoignition  or  hypergolicity  is  therefore  of  major  interest  and  a  significant 
challenge  in  finding  a  suitable  MMH  replacement.  Only  in  the  last  15  years  has  there  been 
much  effort  to  fundamentally  understand  the  hypergolic  chemistry  of  MMH/NTO,  even 
though  the  utility  of  this  combination  has  been  known  since  the  1950s.  The  first  process  in 
finding  a  suitable  hypergol  replacement  has  been  in  identifying  fuels  that  are  environmentally 
friendly,  or  at  least  easier  to  handle,  store  and  transport.  Certain  ionic  liquids  have  been 
shown  to  provide  these  advantages.  The  potential  to  develop  hypergolic  ILs  as  a  suitable 
option  for  space  propulsion  applications  requires  the  availability  of  comprehensive  gas  phase 
thermodynamics  and  chemical  kinetics  databases.  The  prospect  of  hypergolically  igniting 
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these  ILs  also  requires  a  detailed  understanding  of  the  underlying  chemical  pathways  and  the 
reaction  mechanisms  involved  when  ILs  are  treated  with  NTO  or  any  other  suitable  oxidizers 
such  as  inhibited  red  fuming  nitric  acid  (IRFNA)  or  white  fuming  nitric  acid  (WFNA).  The 
literature  in  this  regards  is  very  sparse  for  ILs.  Therefore  the  ignition  chemistry  and 
combustion  chemistry  of  these  ILs  remain  largely  unknown  as  is  their  hypergolicity  with 
NTO,  IRFNA,  WFNA  or  hydrogen  peroxide. 


2.  General  description  of  ionic  liquids  (ILs) 

ILs  are  synthesized  by  reaction  of  a  base  with  a  strong  acid.  ILs  are  therefore  salts  and  consist 
of  negative  and  positive  ions.  It  is  generally  accepted  that  the  anion  (negative  ion)  is  more 
important  in  controlling  the  physical  properties  of  the  IL  and  that  the  nature  of  the  cation 
(positive  ion)  affects  the  stability  and  reactivity  of  the  IL.  Therefore,  it  is  possible  to  adjust 
physical  properties  by  making  changes  to  the  structure  of  either  the  anion  or  the  cation,  or 
both,  to  suit  requirements.  Such  strategies  may  also  be  of  interest  for  designing  energetic 
materials.  Approximately  1018  different  ionic  liquids  are  possible  (Trohalaki  et  al.,  (2005)). 
The  number  of  ILs  synthesized  has  expanded  almost  exponentially  over  the  past  two  decades 
and  it  is  certain  that  new  classes  will  be  developed  in  the  forthcoming  years,  hence  it  is  not 
practical  to  study  experimentally  the  chemical  activity  of  all  of  these  innumerable  ILs. 
Therefore,  theoretical  tools  have  been  developed  for  some  of  the  properties,  and  new  tools 
need  to  be  developed  to  predict  remaining  physicochemical  properties  of  interest  so  as  to 
allow  for  an  efficient  approach  to  screen  and  design  ILs  for  the  desired  application. 
Applications  of  ILs  are  also  numerous:  electro-deposition  of  metals,  fuel  desulphurization, 
propellants,  fuel  cells,  etc.  However,  to  date  only  a  handful  of  the  proposed  applications  have 


5 


come  to  full  fruition.  By  far  the  largest  academic  activity  with  ILs  has  been  with  their 
chemical  synthesis.  Data  on  some  of  the  physiochemical  properties  of  ILs  can  be  found  in  our 
interim  report  #  1 . 


3.  Chemical  nature  of  vapor  above  ILs  and  above  ILs-white  fuming  nitric  acid  (WFNA) 
mixtures 

The  chemical  nature  of  the  vapor  above  an  ionic  liquid  is  difficult  to  determine 
experimentally.  It  is,  however,  important  to  know  how  the  vapor  pressure  forms  over  an  ionic 
liquid  surface.  Four  possibilities  exist  at  equilibrium  above  the  liquid  phase;  an  ion-pair 
(including  ion-pair  clusters),  separated  ions,  neutral  pairs  or  separated  neutrals.  Experimental 
findings  and  theoretical  calculations  are  reported  in  the  literature  on  this  issue.  Ionic  liquids 
can  be  divided  into  two  classes:  protic  ionic  liquids  and  aprotic  ionic  liquids. 

3.1  Protic  ionic  liquids 

According  to  Earle  et  al.  (2006),  an  equilibrium  exists  between  ionic  and  non-ionic  species 
for  protic  ionic  liquids: 

BH+X(£)  B(£)  +  HX(e)  <-+  B(g)  +  HX(g) 

where  (f)  represents  the  liquid  phase  and  (g)  represents  the  gas  phase.  Calculations  made  by 
Schmidt  et  al.  (2005)  on  protic  triazolium-based  energetic  ionic  liquids  show  that  the  neutral 
pair  is  of  lower  energy  than  the  ionic  pair.  However,  separation  of  the  molecules  in  the  pair 
requires  a  significant  amount  of  energy:  about  50  kJ  mol'1.  One  of  the  conclusions  of  this 
paper  is  that  the  ionic  pair  once  formed  will  transfer  one  proton  to  form  a  neutral  pair.  This 
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proton  transfer  is  therefore  fundamental  for  understanding  the  decomposition  of  protic 
triazolium-based  species.  However,  Baranyai  et  al.  (2004)  do  not  observe  the  formation  of 
species  normally  associated  with  the  proton  transfer  process.  Zom  et  al.  (2006)  deduced  from 
their  calculations  on  protic  tetrazolium-based  ionic  liquids  that  when  a  cation  is  paired  with  an 
oxygen  rich  anion,  an  isolated  gas  phase  ion-pair  was  not  generally  found  to  be  stable.  In  fact, 
a  proton  transfers  without  barrier  from  the  cation  to  the  anion  to  form  a  neutral  pair.  This 
result  is  consistent  with  the  experimental  findings  of  Leal  et  al.  (2007),  who  report  that  the  gas 
phase  above  a  protic  ionic  liquid  is  composed  of  separated  neutral  molecules. 

3.2  Aprotic  ionic  liquids 

Calculations  made  by  Emel’yanenko  et  al.  (2007)  on  the  gas  phase  equilibrium: 

ion-pair  <->  separated  ions 

for  l-butyl-3-methylimidazolium  dicyanamide  indicate  that  the  dissociation  of  the  ion-pair  is 
marginal.  Experimental  and  theoretical  works  reported  by  Strasser  et  al.  (2007)  suggest  that 
the  vapor  pressure  of  an  aprotic  ionic  liquid  (1 -ethyl-3 -methylimidazolium 
bis(trifluoromethylsulfonyl)imide)  is  composed  of  isolated  ion  pairs.  The  presence  of  clusters 
of  ion  pairs  also  seems  to  be  negligible.  This  conclusion  is  supported  by  another  experimental 
work  performed  by  Leal  et  al.  (2007),  who  find  that  the  gas  phase  is  composed  of  isolated  ion 
pairs  and  that  no  experimental  evidence  supports  the  formation  of  ion  clusters.  Armstrong  et 
al.  (2007)  performed  the  analysis  of  the  vapor  pressure  of  1 -ethyl-3 -methylimidazolium 
ethylsulfate  by  mass  spectrometry  (MS).  They  deduced  from  their  MS  experiments  that  the 
vapor  of  these  ionic  liquids  consists  of  ion  pairs,  although  they  do  not  observe  directly  these 
ion  pairs  because  the  ionized  ion  pairs  (formed  in  the  ionization  chamber  of  the  mass 
spectrometer  by  electron  bombardment)  are  thought  to  be  unstable  and,  therefore,  undergo 
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fragmentation.  One  other  finding  of  Armstrong  et  al.  (2007)  is  that  clusters  of  ion  pairs  do  not 
seem  to  form  in  the  gas  phase. 

3.3  Chemical  nature  of  gases  above  ILs-WFNA  liquid  mixtures 

One  can  expect  to  observe  in  the  gas  phase,  nitric  acid  (with  small  amounts  of  NO2)  and  other 
neutral  molecules  (neutral  pairs  or  separated  neutrals)  according  to  subsections  3.1  and  3.2. 
Actually,  Chambreau  and  Vaghjiani  et  al.  (2008)  showed  that  ILs  and  WFNA  react 
exothermically  in  the  liquid  phase  to  form  products  in  the  liquid  phase  that  then  elute  into  the 
gas  phase.  This  is  further  discussed  in  section  5. 


4.  Thermochemistry  of  imidazoles,  triazoles  and  tetrazoles 

Imidazole,  triazole,  and  tetrazole  compounds  and  their  derivatives  are  interesting  molecules 
because  of  their  potential  applications  as  explosives  or  propellants.  As  shown  in  section  3, 
these  compounds  are  potentially  interesting  in  the  field  of  ionic  liquids  (ILs)  because 
imidazolium,  triazolium,  and  tetrazolium  cations  are  frequently  encountered  in  ILs.  In  fact, 
even  for  the  simplest  of  these  compounds,  there  is  a  general  lack  of  good  thermodynamics 
data  sets  that  include  the  standard  enthalpy  of  formation  at  298.15  K,  heat  capacity  (at 
constant  pressure,  Cp°  or  at  constant  volume,  Cv°)  for  a  wide  range  of  temperatures,  and  the 
standard  absolute  entropy,  S°.  Such  data  is  needed  to  carry  out  accurate  kinetics  modeling 
studies. 


8 


4. 1  Estimation  of  thermodynamic  data 


Ab-initio  methods  have  been  shown  to  be  reliable  tools  for  estimating  the  thermochemistry  of 
many  compounds.  Three  methods  of  estimating  the  gas  phase  standard  enthalpy  of  formation 
at  298.15  K  have  been  proposed  by  the  author  of  this  report  using  quantum  chemistry  tools. 
Two  of  these  methods  can  be  applied  to  imidazole,  triazole  and  tetrazole  compounds.  One  of 
these  methods  is  devoted  to  general  organic  compounds  and  the  other  to  energetic  materials 
with  nitro  and  nitrate  groups.  The  gas  phase  standard  enthalpy  of  formation  of  molecule  j  at 
298.15  K  can  be  determined  from  the  following  equation: 

(  \ 


AfH  298.15 K(g) 


=  627.51  x 


EJ 


ZPEj  +  thermal  corrections  +  ^  ai  c* 


V  i  J 

where  at  is  the  number  of  atoms  i  in  molecule  j  and  c*  the  atomic  correction  for  atom  i.  Ej 

and  ZPEj  denote,  respectively,  the  absolute  electronic  energy  and  zero-point  energy, 
calculated  using  either  Gaussian  98W  or  GAUSSIAN  03  suite  of  programs.  The  units  are 
Hartree  molecule"1  for  Ej,  ZPEj  and  thermal  corrections,  and  Hartree  atom'1  for  c* ,  whereas 


AfH029s.\5K(g)  is  in  kcal  mol"1.  Three  atomic  corrections  (C,  H  and  N)  are  considered  here 

for  the  method  established  for  organic  compounds  and  five  atomic  corrections  (C,  H,  N,  Cd, 
Nd)  are  considered  for  the  method  established  for  energetic  compounds  with  one  more  atomic 
correction  for  hypervalent  nitrogen  (noted  Nhyp)  as  encountered  in  nitroimidazoles, 
nitrotriazoles  and  nitrotetrazoles.  Cp°,  Cv°,  and  S°  are  calculated  using  statistical 
thermodynamic  equations  once  the  vibration  frequencies  and  inertia  moments  have  been 
estimated  using  the  Gaussian  software  packages. 
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4.2  Results 


Table  1  compares  the  experimental  data  and  the  calculated  data  with  both  methods,  and  Table 
2  reports  calculated  gas  phase  standard  enthalpy  of  formation  at  298.15  K  with  both  methods 
for  compounds  for  which  no  experimental  data  is  available. 


Table  1.  Comparison  between  experimental  gas  phase  standard  enthalpy  of  formation  at 
298.15  K.  and  calculated  gas  phase  standard  enthalpy  of  formation  at  298.15  K  for 
imidazoles,  triazoles  and  tetrazoles. _ 


Compound 

CAS  Registry 
Number 

Experimental 

A  fH° 298A5K  (g) 

(kcal  mol"1) 

Calculated 

AfH029$.\5K(g)  with 

the  method  for 
organic  compounds 
(kcal  mol"1) 

Calculated 

AfH029SA5K(g)  with 

the  method  for 
energetic  compounds 
(kcal  mol"1) 

1H- 1,2,4-triazole 
288-88-0 

N - T\ 

O 

N 

H 

46.1±0.2 

45.05 

40.7 

43.4 

2-ethyl- 1  H-imidazole 
1072-62-4 

/-O' 

H3C  NT 

H 

16.3±0.31 

13.5 

15.6 

lH-tetrazole 

288-94-8 

O 

N - N 

76.6T0.7 

71.7 

74.0 

1H- imidazole 
288-32-4 

\J 

H 

33.29i0.45 

30.6T1.8 

29.1 

32.3 
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2-methyl- 1H- 
imidazole 
693-98-1 

H 

21.5±0.26 

18.5 

21.1 

lH-l,2,3-triazole 

288-36-8 

O 

^ - N 

59 

57 

59.7 

lH-pyrazole 

288-13-1 

C> 

42.9 

43.3 

39.2 

42.4 

Pyrrole 

109-97-" 

H 

N 

i 

25.9T0.1 

34.23 

26.4 

30.0 

Table  1  shows  unambiguously  that  the  method  established  for  energetic  compounds  performs 
better  than  the  method  established  for  general  organic  compounds.  This  conclusion  is 
interesting  because  this  method  has  been  established  for  nitro  and  nitrate  compounds  and 
therefore  it  can  be  expected  to  provide  reasonably  accurate  results  for  nitroimidazoles, 
nitrotriazoles  and  nitrotetrazoles  reported  in  Table  2.  Calculations  are  computationally 
expensive  and  all  of  the  results  are  not  available  at  this  time  for  most  of  the  compounds 
reported  in  Table  2.  Because  of  the  importance  of  these  compounds  in  the  field  of  energetic 
materials,  a  complete  data  set  of  these  results  will  be  published  in  a  journal  article  at  a  later 
date. 
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Table  2,  Calculated  gas  phase  standard  enthalpy  of  formation  at  298.15  K  for  imidazole. 


triazole  and  tetrazole  compounds  for  which  experimental  data  are  not  avai 

able. 

Compound 

CAS  Registry  Number 

Experimental 

A  y//°298.15X(g) 

(kcal  mol"1) 

Calculated 

A^//°298.15X(g) 

with  the  method 
for  organic 
compounds  (kcal 
mol'1) 

Calculated 

A y’/i/0 298.15  X  (g) 

with  the  method 
for  energetic 
compounds  (kcal 
mol'1) 

1  -methyl- 1  H-imida2 
616-47-7 

H3C — 

x>le 

N 

26.2 

28.9 

1 ,2-dimethyl- 1 H- imidazole 
1739-84-0 

ch3 

Or- 

^ - N 

16.3 

18.4 

4-methylimidazol 

822-36-6 

h3c - <(  ] 

e 

H 

19.6 

22.2 

l-Methyl-5- 

nitroimidazoh 

3034-42-2 

yt 

4-methyl-5- 

nitroimidazob 

14003-66-8 

A  , 

/ 

1  1 

3  -nitro- 1 H- 1 ,2,4-triazole 
24807-55-4 

■Or- 

1 ,5-diamino- 1 ,2,3 ,4-tetrazole 

90.6 

99.8 

4-amino- 1 ,2,4-triazole 

73.7 

79.8 

5.  Thermochemistry  of  compounds  formed  during  ILs-WFNA  reaction 

Chambreau  and  Vaghjiani  et  al.  (2008)  recently  presented  experimental  results  on  a  class  of 
ILs  which  exhibit  hypergolicity  at  room  temperature  with  WFNA.  They  presented 
mechanistic  information  on  reactions  in  the  liquid  phase  to  explain  the  formation  of  species 
observed  by  Fourier  transform  infrared  (FTIR)  spectroscopy  in  the  gas  phase  above  the 
reacting  ILs  and  WFNA  mixtures.  For  1 -butyl-3 -methyl  imidazolium  dicyanamide  reacting 
with  WFNA  they  observed  in  the  gas  phase;  HNO3,  CC^,  N2O  and  HNCO.  Also  present  was 
water  and  NO2  (WFNA  contains  some  amount  of  NO2).  HNCO,  N2O  and  CO2  were 


presumed  to  form  from  dinitrobiuret,  NO2NHCONHCONHNO2  (abbreviated  DNB), 


which  decomposes  thermally  to  give  NO2NHCONCO  (or  intermediate  #  1,  see  below). 


Sketch  2:  molecular  structure  of  intermediate  #1 


and  NH2NO2  (or  intermediate  #  2,  see  below). 
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Sketch  3:  molecular  structure  of  intermediate  #2 


These  two  intermediates  in  the  liquid  phase,  respectively,  give  HNCO  +  N20  +  C02  and  N20 
+  H20.  These  compounds  together  with  HNO3  and  N02  from  the  WFNA  will  lead  to  ignition 
or  not  depending  on  the  temperature  and  the  concentration  of  the  species  reached  in  the  gas 
phase.  Since  DNB,  intermediate  #1  and  intermediate  #2  are  nitro  compounds,  their 
thermochemistry  can  be  calculated  by  using  the  theoretical  method  established  for  nitro  and 
nitrate  compounds.  However,  these  three  compounds  are  not  observed  in  the  gas  phase.  It  is 
therefore  probably  more  logical  to  tentatively  explain  ignition  by  considering  only  those 
species  that  are  actually  observed  in  the  gas  phase  namely;  HNCO,  N20,  N02,  HNO3,  H20 
and  C02. 

5.1  Gas  phase  thermochemistry  of  DNB  and  of  the  intermediates 

The  gas  phase  thermochemistry  of  DNB  can  be  calculated  by  using  the  method  established  for 
nitro  compounds.  The  validation  of  this  method  was  presented  by  Osmont  et  al.  (2007).  This 
procedure  gives  a  value  of  -47.9  kcal  mol"1  for  the  gas  phase  standard  enthalpy  of  formation  at 
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298. 15K.  Vibration  frequencies  and  inertia  moments  have  been  computed  to  allow  for  the 
calculation  of  Cp°  and  S°  over  a  wide  temperature  range.  For  intermediate  #  2,  or  nitramide, 
the  experimental  data  available  are  somewhat  confusing.  Burcat  and  Ruscic  (2009)  in  the 
latest  version  of  their  thermodynamic  database  report  two  quite  different  values  for  the  gas 
phase  standard  enthalpy  of  formation  at  298.15  K:  -6.2±2.4  kcal  mol"1  and  -0.7  kcal  mol"1. 
Our  method  applied  to  NH2NO2  leads  to  a  value  of  +1.4±3.0  kcal  mol"1,  which  is  more 
consistent  with  their  latter  number.  Vibration  frequencies  and  inertia  moments  have  been 
computed  to  allow  for  the  calculation  of  Cp°  and  S°  over  a  wide  temperature  range.  For 
intermediate  #  1,  a  value  of  -44.5  kcal  mol"1  is  derived  by  using  the  same  method  as  that  for 
DNB  and  nitramide.  As  with  the  other  compounds,  vibration  frequencies  and  inertia  moments 
have  been  computed  for  intermediate  #  1  to  allow  for  the  calculation  of  Cp°  and  S°  over  a 
wide  temperature  range. 

5.2  Thermochemistry  of  HNCCX  HNCK  NO?,  N?Q,  H?Q,  and  CO? 

Species  present  as  reactants  in  the  gas  phase  above  reacting  liquids  are  HNCO  (isocyanic 
acid),  nitric  acid  (HNO3  or  HONO2X  N02,  N2O,  H20  and  CO2.  The  thermochemistry  of  these 
small  species  are  generally  known  and  are  reported  in  the  CHEMKIN  format  by  Burcat  and 
Ruscic  (2009)  in  their  latest  version  of  their  thermodynamics  database  and  are  also  reported 
below  in  the  CHEMKIN  format. 

HNCO  Isocyanic  AciA  5/05H  l.N  l.C  1.0  l.G  200.000  6000.000  B  43.02478  1 

5 . 30045051E+00  4 . 02250821E-03-1 . 40962280E-06  2 . 23855342E-10-1 . 32499966E-14  2 

-1 . 61995274E+04-3.11770684E+00  2 . 24009031E+00  1 . 45600497E-02-1 . 54352330E-05  3 

8.55535028E-09-1.79631611E-12-1 . 54589951E+04  1 . 21 66377 5E+01-1 . 42 64274 0E+04  4 

HN03  T  8/03H  l.N  1.0  3.  0.G  200.000  6000.000  B  63.01288  1 
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8 . 03098942E+00  4 . 46958589E-03-1 . 72459491E-06  2 . 91556153E-10-1 . 80102702E-14  2 

-1 . 93138183E+04-1 . 62616537E+01  1 . 69329154E+00  1 . 90167702E-02-8 . 2517 6697E-06  3 

-6.06113827E-09  4 . 65236978E-12-1 . 74198909E+04  1 . 71839838E+01-1 . 61524852E+04  4 

N02  L  7/88N  10  2  0  OG  200.000  6000.000  A  46.00554  1 

0 . 48847540E+01  0 . 21723955E-02-0 . 82806909E-06  0 . 1574 751 OE-O 9-0 . 1051 08 95E-13  2 

0 . 23164982E+04-0 . 11741695E+00  0 . 39440312E+01-0 . 15854290E-02  0 . 16657812E-04  3 

-0 .20475426E-07  0 . 78350564E-11  0 . 28966180E+04  0 . 6311 991 9E+01  0 . 41124701E+04  4 

N20  L  7/88N  20  1  0  OG  200.000  6000.000  A  44.01288  1 

0 .48230729E+01  0 . 26270251E-02-0 . 95850872E-06  0 . 160 007 12E- 09-0 . 97752302E-14  2 

0. 80734 047E+04-0. 22 0172 08E+01  0 . 22571502E+01  0 . 11304728E-01-0 . 1367 131 9E- 04  3 

0 . 96819803E-08-0 . 29307182E-11  0 . 87417746E+04  0 . 10757992E+02  0 . 98141682E+04  4 

H20  L  5/89H  20  1  0  OG  200.000  6000.000  A  18.01528  1 

0 .26770389E+01  0 . 29731816E-02-0 . 77376889E-06  0 . 9443351 4E-10-0 . 4268 9991E-14  2 

-0 .29885894E+05  0 . 68825500E+01  0 . 41986352E+01-0 . 20364017E-02  0 . 65203416E-05  3 

-0.54879269E-08  0 . 1771968 0E-11-0 . 302 9372 6E+05-0 . 84900901E+00-0 . 29084817E+05  4 

C02  L  7/88C  10  2  0  OG  200.000  6000.000  A  44.00980  1 

0 .46365111E+01  0 . 27414569E-02-0 . 99589759E-06  0 . 16038 66 6E-0 9-0 . 91619857E-14  2 

-0 . 49024904E+05-0 . 19348955E+01  0 . 23568 130E+01  0 . 8984129 9E-02-0 . 71220 632E-05  3 

0 .24573008E-08-0 . 14288548E-12-0 .48371971E+05  0 . 9900 9035E+01-0 . 4732 8105E+05  4 


6.  Calculations  at  thermodynamic  equilibrium 


For  non-hypergolic  systems,  i.e.  systems  for  which  an  external  ignition  device  is  required, 
chemical  equilibrium  calculations  are  often  used  to  estimate  flammability  limits.  From  these 
equilibrium  calculations,  adiabatic  flame  temperatures  are  considered  as  a  criterion  to  estimate 
lower  and  upper  flammability  limits,  at  a  given  total  pressure.  It  has  been  shown  that  the 
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adiabatic  flame  temperature  at  the  lower  and  upper  flammability  limits  at  a  given  pressure 
remains  approximately  constant  in  hydrocarbon/air  mixtures  and  the  reported  experimental 
limits  at  atmospheric  pressure  correspond  generally  to  adiabatic  flame  temperatures  in  the 
1000-1500  K  range.  A  criterion  of  1200  K  is  often  found  convenient.  It  is  therefore  tenable 
that  the  flammability  diagrams  of  hypergolic  mixtures  can  be  rationalized  with  such  a 
thermochemical  criterion.  All  these  points  were  discussed  in  interim  report  #2.  A  criterion 
of  about  2000  K  and  above  seems  to  be  appropriate  to  discuss  the  occurrence  of  hypergolic 
ignition  in  a  system.  This  simple  thermochemical  criterion  is  actually  not  sufficient  and  needs 
to  be  complemented  by  a  detailed  chemical  kinetics  model  to  prove  that  a  system  can  be 
hypergolic  under  specified  conditions.  Concerning  ionic  liquids  and  WFNA,  ignition  has 
been  observed  in  a  gas  phase  mixture  of  HNCO,  N2O,  NO2,  HNO3,  H2O  and  CO2  at 
unspecified  conditions  of  temperature  and  pressure.  The  first  task  is  to  ascertain  whether  a 
mixture  containing  these  species  has  the  ability  to  sustain  a  flame.  In  the  framework  of 
interim  report  #  2,  this  task  was  developed  in  general  terms  and  then  specifically  applied  to 
the  MMH/NTO  hypergolic  system.  Such  an  analysis  can  also  be  performed  here,  but  as  in  the 
case  of  MMH/NTO,  the  exact  composition  of  the  mixture  in  which  ignition  occurs  is  more 
than  often  not  well  characterized.  Future  studies  at  AFRL,  Edwards  AFB  to  determine  the 
composition  of  such  mixtures  would  be  very  valuable  in  the  validation  procedures  discussed 
below. 

Since  the  quantitative  composition  of  the  gas  phase  above  the  liquid  is  not  known,  and  even  if 
it  can  be  considered  qualitatively  as  known,  the  best  approach  to  use  here  is  to  examine  many 
mixtures  at  different  initial  conditions  of  pressure  and  temperature  to  determine  if  such 
mixtures  are  able  to  sustain  a  flame.  The  answer  is  necessarily  yes,  because  ignition  is 
experimentally  observed  (Chambreau  and  Vaghjiani  et  al.  (2008)).  It  remains  possible, 
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however,  that  an  important  species  needed  to  quantitatively  explain  the  observed  ignition 


delay  is  presently  unidentified  and  its  nature  must  be  identified  to  amend  the  initial  detailed 
chemical  kinetics  model.  For  instance,  species  such  as  O2  cannot  be  observed  by  using  the 
FTIR  technique. 


6.1  Pure  HNCO 


A  determination  of  whether  pure  HNCO  is  able  to  sustain  a  flame  is  an  important  exercise. 
The  gas  phase  standard  enthalpy  of  formation  at  298.15  K  of  HNCO  is  about  -28  kcal  mol"1. 
This  would  generally  indicate  that  the  decomposition  of  such  a  molecule  would  be 
endothermic.  This  is,  however,  not  true  for  some  compounds,  as  is  the  case  with  certain 
organic  peroxides  and  perchloric  acid  (HCIO4).  Equilibrium  calculations  can  shed  some  light 
in  this  area.  For  pure  HNCO,  adiabatic  equilibrium  calculations  at  constant  pressure  indicate 
that  its  thermal  decomposition  is  exothermic  as  shown  below: 


INITIAL  STATE: 

EQUILIBRIUM  STATE 

P  (atm) 

1 . 0000E  +  00 

1 . 0000E+00 

T  (K) 

3 . 8000E+02 

8 . 8078E+02 

Mole  Fractions 

HNCO 

1 . 0000E+00 

8.1567E-08 

CO 

0.000  0E  +  0  0 

4.17  68E-01 

C02 

0.000  0E  +  0  0 

8 . 2333E-02 

H20 

0.0000E+00 

8 . 2704E-03 

CH4 

0.000  0E  +  0  0 

9 . 0590E-02 

H2 

0.000  0E  +  0  0 

1 . 0580E-01 

N2 

0.000  0E  +  0  0 

2 . 9529E-01 
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At  equilibrium,  a  temperature  of  about  880  K  is  predicted.  Since  this  is  lower  than  the  2000  K 
adiabatic  flame  temperature  criterion,  this  probably  means  that  pure  HNCO  is  not  able  to 
sustain  a  decomposition  flame.  However,  even  if  the  initial  chosen  state  is  arbitrary,  this 
calculation  shows  that  the  thermal  decomposition  of  HNCO  must  be  considered  in  the 
detailed  chemical  kinetics  model.  The  associated  exothermicity  can  possibly  be  the  cause  of 
ignition  or,  at  least,  an  important  parameter  to  consider.  This  can  only  be  examined  once  a 
detailed  chemical  kinetics  model  is  built. 


6.2  Mixtures  of  HNCO  and  HNO3 


Here  we  examine  the  possibility  that  ignition  is  specifically  due  to  HNCO-HNO3  mixtures.  A 
50/50  mixture  is  chosen  as  a  starting  example. 


INITIAL  STATE: 

EQUILIBRIUM  STATE 

P  (atm) 

1 . 0000E+00 

1 . 0000E+00 

T  (K) 

3 . 8000E+02 

2 . 5090E+03 

Mole  Fractions 

HNCO 

5 . 0000E-01 

4 . 9302E-10 

CO 

0.00  0  0E  +  0  0 

2 . 5932E-02 

O 

0.0000E+00 

5 . 6990E-03 

H 

0.0000E+00 

1.7171E-03 

C02 

0.00  0  0E  +  0  0 

2 . 5258E-0 1 

NO 

0.0000E+00 

1 . 1339E-02 

H20 

0.0000E+00 

2 . 6137E-0 1 

HN03 

5 . 0000E-01 

2 . 9272E-11 

OH 

0.0000E+00 

2 . 4029E-02 

H2 

0.0000E+00 

4 . 2509E-03 
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N2 

0.00  0  0E  +  0  0 

2 . 7283E-01 

02 

0.00  0  0E  +  0  0 

1.4022E-01 

H02 

0.0000E+00 

1 . 9371E-05 

The  adiabatic  flame  temperature  at  constant  pressure  is  about  2500  K.  This  indicates  that 
such  a  system  is  able  to  sustain  a  flame. 


6.3  Mixtures  of  HNCO,  CO?.  H9O.  andN?Q 


By  considering  the  equations  given  by  Chambreau  and  Vaghjiani  et  al.  (2008)  for  one  mole  of 
DNB  decomposed  in  the  liquid  phase,  there  will  be  an  equal  amount  of  HNCO,  CO2,  H2O  and 
twice  this  amount  of  N2O  in  the  gas  phase.  Such  a  mixture  can  probably  lead  to  hypergolic 
ignition  as  shown  below  because  an  adiabatic  temperature  of  about  2300  K  is  computed  in 
which  C02,  H20,  N2  and  O2  represent  about  97  mol%  of  the  gases  at  thermodynamic 
equilibrium. 


P  (atm) 
T  (K) 


INITIAL  STATE: 
1 . 0000E  +  00 
3 . 8000E+02 


EQUILIBRIUM  STATE: 
1 . 0000E  +  00 
2 . 2  833E  +  03 


Mole  Fractions 


HNCO 

2 . 0000E-01 

2 . 4181E-10 

N20 

4 . 0000E-01 

1 . 9379E-07 

CO 

0.0000E+00 

1 . 5069E-02 

O 

0.000  0E  +  0  0 

9 . 4 1 94E-04 

H 

0.000  0E  +  0  0 

4.047  0E-04 

C02 

2 . 0000E-01 

3 . 0131E-01 

NO 

0.000  0E  +  0  0 

4 . 9068E-03 
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H20 

2 . 0000E-01 

2.3109E-01 

OH 

0.00  0  0E  +  0  0 

7.87  83E-03 

H2 

0.0000E+00 

2.0522E-03 

N2 

0.00  0  0E  +  0  0 

3. 9302E-01 

02 

0.0000E+00 

4 . 3310E-02 

N2O  is  known  to  be  an  oxidant  as  well  as  a  monergol,  i.e.  its  decomposition  is  sufficiently 
exothermic  that  it  is  able  to  sustain  a  decomposition  flame.  Therefore,  the  thermal 
decomposition  of  N2O  must  be  included  in  the  detailed  chemical  kinetics  model.  It  can  be 
supposed  that  for  the  present  case,  thermal  decomposition  of  N2O  leads  to  the  production  of  O 
atoms,  through  the  spin-forbidden  process,  N20  —>  N2  +  O,  and  that  ground-state  O  atoms 
then  further  react  with  HNCO.  The  former  process  is  known  to  be  the  initiation  step  during 
thermal  decomposition  of  pure  N2O. 


6.4  Mixtures  of  HNCO,  CO 2.  H9O.  N9O.  HNCh.  and  NO 2 


The  amount  of  HNO3  in  the  gas  phase  will  be  related  to  the  concentration  of  HNO3  in  the 
liquid  solution  of  IL  +  WFNA  at  a  given  temperature.  Raoult’s  law  can  be  used  for  ideal 
solutions  or  activity  coefficients  can  be  used  for  non-ideal  solutions  to  predict  the 
composition.  It  should  be  noted  that  a  mixture  of  WFNA  and  water  was  also  used  in  some  of 
the  experiments  performed  by  Chambreau  and  Vaghjiani  et  al.  (2008).  Moreover,  small 
amounts  of  NO2  are  also  present  in  WFNA.  Presence  of  HNO3  and  NO2  in  the  mixture  allows 
for  two  oxidants  to  be  considered. 

INITIAL  STATE:  EQUILIBRIUM  STATE: 

P  (atm)  1 . 0000E+00  1.0000E+00 

T  (K)  3 . 8000E+02  1.9195E+03 
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Mole  Fractions 


HNCO 

1 . 6393E-01 

1 . 1061E-12 

CO 

0.000  0E  +  0  0 

3 . 5231E-04 

N02 

1 . 6393E-02 

7 . 1443E-06 

N20 

3.27  87E-01 

1 . 4522E-07 

0 

0.000  0E  +  0  0 

1.5090E-04 

C02 

1 . 6393E-01 

2 . 3066E-01 

NO 

0.000  0E  +  0  0 

3 . 8174E-03 

H20 

1 . 6393E-01 

2 . 2  98  OE-0 1 

HN03 

1 . 6393E-01 

4.207  0E-1 1 

OH 

0.000  0E  +  0  0 

2 . 2  64  0E-03 

N2 

0.000  0E  +  0  0 

3.5039E-01 

02 

0.000  0E  +  0  0 

1 . 8247E-01 

In  this  case,  the  adiabatic  flame  temperature  at  constant  pressure  is  about  2000  K  and  CCE, 
H2O,  N2  and  O2  constitute  about  99.3  mol%  of  the  gases  at  thermodynamic  equilibrium.  A 
mixture  with  this  composition  is  probably  able  to  sustain  a  flame. 


7.  Detailed  chemical  kinetics  model 


7.1  Introduction 


It  seems  clear,  from  the  analyses  presented  in  section  6,  that  the  species  identified  in  the  gas 
phase  namely;  HNCO,  HNO3,  N2O,  NO2,  CO2,  and  H2O  are  able  to  sustain  a  flame. 
However,  thermodynamics  calculations  at  equilibrium  alone  do  not  explain  by  which 
chemical  pathway (s)  ignition  occurs.  Elucidation  of  these  pathways  is  relevant  to  the 
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development  of  a  detailed  chemical  kinetics  model,  and  such  a  model  is  needed  to  understand 
the  chemical  reactivity  of  a  system  when  time  is  a  parameter  of  the  problem.  Data  of  interest 
which  can  be  estimated  with  a  detailed  chemical  kinetics  model  are  ignition  times  and  species 
time  profiles  during  the  ignition  delay.  Species  profiles  may  also  be  determined  when  there  is 
no  ignition  in  the  system.  In  fact,  depending  on  several  other  parameters  such  as  composition, 
temperature,  and  pressure,  there  will  be  a  finite  time  to  reach  equilibrium  from  a  given  initial 
state  of  the  system  under  consideration. 

7.2  Building  of  a  detailed  chemical  kinetics  model 

It  is  beyond  the  scope  of  this  work  to  present  all  of  the  features  and  considerations  of  the  field 
of  chemical  kinetics.  Some  of  the  important  facts  are  discussed  below.  A  detailed  chemical 
kinetics  model  is  composed  of  many  reactions  depending  on  the  number  of  reactants  and  the 
number  of  atoms  in  the  molecules.  The  chemical  nature  of  the  atoms  involved  is  also 
important.  All  these  reactions  are  considered  to  be  elementary.  They  can  be  unimolecular, 
bimolecular  or  termolecular.  Unimolecular  and  termolecular  reactions  are  pressure-dependent 
and  there  are  several  mathematical  methods  to  treat  them.  A  detailed  chemical  kinetics  model 
will  have  initiation  reactions,  chain  branching  reactions,  chain  propagation  reactions  and  chain 
termination  reactions.  From  the  results  of  the  previous  section,  it  can  be  concluded  that  the 
thermal  decomposition  of  HNCO,  HNO3  and  N20  must  be  written  in  detail  to  account  for  the 
observed  exothermicity,  and  also  because  of  the  subsequent  production  of  additional  reactive 
species  in  these  processes.  The  atoms  and  radicals  formed  during  the  three  corresponding 
initiation  reactions  will  interact  with  the  remaining  HNCO,  HNO3,  and  N20,  and  also  with  the 
any  other  species  present  in  the  system  such  as  H20,  C02  and  N02.  The  radical  species,  N02 
will  also  be  available  for  reactions  in  the  system.  Once  all  of  the  feasible  reactions  are 
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written,  at  least  as  best  as  possible,  a  rate  constant,  k,  needs  to  be  specified  for  each  of  the 
reaction  steps.  Depending  on  the  availability  of  thermodynamics/kinetics  databases  a  value  of 
k  may  be  uniquely  assignable,  k  is  generally  considered  to  follow  the  Arrhenius  equation: 

-E 

k  =  Ae«r 

where  A  represents  the  frequency  (or  pre-exponential)  factor,  E  the  activation  energy,  R  the 
perfect  gas  constant  and  T  the  temperature.  However,  the  A  factor  also  has  a  temperature 
dependence,  which  generally  cannot  be  neglected  for  accurate  combustion  predictions  and 
therefore,  a  modified- Arrhenius  equation  is  often  needed: 

-E 

k-  A  Tn  e™ 

The  following  pages  report  our  initial  detailed  chemical  kinetics  model  relevant  to  the  present 
system.  The  model  consists  of  61  species  and  377  reactions.  Rate  constants  are  from 
available  compilations  and  in  some  cases  have  been  estimated  by  analogy  with  similar 
reactions.  Units  used  are;  E(cal  mol'1),  T(K),  R(cal  mol'1  K'1),  and  A(cm3  mol'1  s'1). 


c 

P  H 
H  A 
A  R 

SPECIES  S  G  MOLECULAR  TEMPERATURE  ELEMENT  COUNT 


CONSIDERED 

E 

E 

WEIGHT 

LOW 

HIGH 

C 

H 

N 

0 

1  . 

HNCO 

G 

0 

43.02522 

200.0 

6000.0 

1 

1 

1 

1 

2  . 

NCO 

G 

0 

42.01725 

300.0 

4000.0 

1 

0 

1 

1 

3. 

N 

G 

0 

14.00670 

200.0 

6000.0 

0 

0 

1 

0 

4  . 

CO 

G 

0 

28.01055 

200.0 

6000.0 

1 

0 

0 

1 

5. 

NO  2 

G 

0 

46.00550 

200.0 

6000.0 

0 

0 

1 

2 

6. 

N20 

G 

0 

44.01280 

200.0 

6000.0 

0 

0 

2 

1 

7  . 

0 

G 

0 

15.99940 

200.0 

6000.0 

0 

0 

0 

1 

8. 

C 

G 

0 

12.01115 

200.0 

6000.0 

1 

0 

0 

0 

9. 

H 

G 

0 

1 .00797 

200.0 

6000.0 

0 

1 

0 

0 

10. 

C02 

G 

0 

44.00995 

200.0 

6000.0 

1 

0 

0 

2 

11. 

NO 

G 

0 

30.00610 

200.0 

6000.0 

0 

0 

1 

1 

12  . 

H20 

G 

0 

18.01534 

200.0 

6000.0 

0 

2 

0 

1 

13. 

NH2N02 

G 

0 

62.02814 

200.0 

6000.0 

0 

2 

2 

2 

14  . 

HN03 

G 

0 

63.01287 

200.0 

6000.0 

0 

1 

1 

3 

15. 

HONO 

G 

0 

47.01347 

200.0 

6000.0 

0 

1 

1 

2 

16. 

NH 

G 

0 

15.01467 

200.0 

6000.0 

0 

1 

1 

0 

17  . 

HNO 

G 

0 

31 .01407 

200.0 

6000.0 

0 

1 

1 

1 

18. 

NH2 

G 

0 

16.02264 

200.0 

3000.0 

0 

2 

1 

0 

19. 

OH 

G 

0 

17.00737 

200.0 

6000.0 

0 

1 

0 

1 
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20. 

CH3 

G 

0 

15.03506 

300.0 

5000.0 

1 

3 

0 

0 

21. 

CH4 

G 

0 

16.04303 

200.0 

6000.0 

1 

4 

0 

0 

22. 

H2 

G 

0 

2.01594 

300.0 

5000.0 

0 

2 

0 

0 

23. 

NH3 

G 

0 

17.03061 

300.0 

5000.0 

0 

3 

1 

0 

24  . 

N2 

G 

0 

28.01340 

300.0 

5000.0 

0 

0 

2 

0 

25. 

CN 

G 

0 

26.01785 

300.0 

5000.0 

1 

0 

1 

0 

26. 

02 

G 

0 

31.99880 

300.0 

5000.0 

0 

0 

0 

2 

27  . 

NO  3 

G 

0 

62.00490 

300.0 

5000.0 

0 

0 

1 

3 

28. 

H02 

G 

0 

33.00677 

300.0 

5000.0 

0 

1 

0 

2 

29. 

H202 

G 

0 

34.01474 

300.0 

5000.0 

0 

2 

0 

2 

30. 

N204 

G 

0 

92.01100 

300.0 

5000.0 

0 

0 

2 

4 

31. 

CH30 

G 

0 

31.03446 

300.0 

3000.0 

1 

3 

0 

1 

32. 

CH 

G 

0 

13.01912 

300.0 

5000.0 

1 

1 

0 

0 

33. 

HCO 

G 

0 

29.01852 

300.0 

5000.0 

1 

1 

0 

1 

34. 

HOCO 

G 

0 

45.01792 

300.0 

4000.0 

1 

1 

0 

2 

35. 

CHCO 

G 

0 

41.02967 

300.0 

4000.0 

2 

1 

0 

1 

36. 

CH2S 

G 

0 

14.02709 

300.0 

4000.0 

1 

2 

0 

0 

37. 

CH2T 

G 

0 

14.02709 

200.0 

6000.0 

1 

2 

0 

0 

38. 

CH20 

G 

0 

30.02649 

300.0 

5000.0 

1 

2 

0 

1 

39. 

CH30H 

G 

0 

32.04243 

300.0 

5000.0 

1 

4 

0 

1 

40. 

CH302 

G 

0 

47.03386 

300.0 

5000.0 

1 

3 

0 

2 

41. 

CH20H 

G 

0 

31.03446 

250.0 

4000.0 

1 

3 

0 

1 

42  . 

CH3CHO 

G 

0 

44.05358 

300.0 

5000.0 

2 

4 

0 

1 

43. 

C2H 

G 

0 

25.03027 

300.0 

4000.0 

2 

1 

0 

0 

44  . 

C2H2 

G 

0 

26.03824 

300.0 

5000.0 

2 

2 

0 

0 

45. 

N2H2 

G 

0 

30.02934 

300.0 

5000.0 

0 

2 

2 

0 

46. 

HCN 

G 

0 

27.02582 

300.0 

4000.0 

1 

1 

1 

0 

47  . 

NNH 

G 

0 

29.02137 

250.0 

4000.0 

0 

1 

2 

0 

48. 

C2H4 

G 

0 

28.05418 

300.0 

5000.0 

2 

4 

0 

0 

49. 

C2H5 

G 

0 

29.06215 

300.0 

5000.0 

2 

5 

0 

0 

50. 

CH2CO 

G 

0 

42.03764 

300.0 

5000.0 

2 

2 

0 

1 

51. 

CH302H 

G 

0 

48.04183 

200.0 

6000.0 

1 

4 

0 

2 

52  . 

C2H3 

G 

0 

27.04621 

300.0 

5000.0 

2 

3 

0 

0 

53. 

CH3CO 

G 

0 

43.04561 

300.0 

5000.0 

2 

3 

0 

1 

54  . 

HCCOH 

G 

0 

42.03764 

300.0 

4000.0 

2 

2 

0 

1 

55. 

CH3NH 

G 

0 

30.04973 

300.0 

5000.0 

1 

4 

1 

0 

56. 

CH2NH 

G 

0 

29.04176 

200.0 

6000.0 

1 

3 

1 

0 

57  . 

C2H6 

G 

0 

30.07012 

300.0 

4000.0 

2 

6 

0 

0 

58. 

H2CN 

G 

0 

28.03379 

300.0 

4000.0 

1 

2 

1 

0 

59. 

HOCN 

G 

0 

43.02522 

300.0 

4000.0 

1 

1 

1 

1 

60. 

C2N2 

G 

0 

52.03570 

300.0 

5000.0 

2 

0 

2 

0 

61. 

HNNO 

G 

0 

45.02077 

300.0 

4000.0 

0 

1 

2 

1 

(k  =  A 

T* *b  exp  ( 

-E/RT) ) 

REACTIONS  CONSIDERED 

A 

b 

E 

1. 

HNCO+M=NH+CO+M 

1.30E+16 

o 

o 

84368.0 

2  . 

HNCO+H=NH2+CO 

3.60E+04 

2.5 

2344.0 

3. 

HNCO+0=HNO+CO 

1.70E+06 

2.1 

o 

o 

4  . 

HNC0+0=NH+C02 

1.70E+06 

2.1 

o 

o 

5. 

HNC0+0H=NH2+C02 

6.30E+10 

i 

o 

i — 1 

11643.0 

6. 

HNCO+OH=NCO+H20 

Declared  duplicate  reaction... 

5.20E+10 

o 

o 

17565.0 

7  . 

HNCO+H=NCO+H2 

1.80E+05 

2.4 

9915.0 

8. 

HNCO+0=NCO+OH 

3.10E+06 

1.9 

6457.0 

9. 

HNCO+OH=NCO+H20 

Declared  duplicate  reaction... 

3.60E+07 

1.5 

3596.0 

10. 

HNCO+CH3=NCO+CH4 

1.00E+12 

o 

o 

9935.0 

11. 

HNCO+NH2=NCO+NH3 

1.00E+12 

o 

o 

8941.0 

12  . 

HNCO+NH=NCO+NH2 

1.00E+12 

o 

o 

8941.0 

13. 

HNCO+N02=NCO+HONO 

1.00E+12 

o 

o 

8941.0 

14  . 

NCO+N02=CO+NO+NO 

2.10E+11 

o 

o 

-874.0 

15. 

NC0+N02=C02+N20 

2.30E+12 

o 

o 

-874.0 

16. 

NCO+OH=H+CO+NO 

8.30E+12 

1 

o 

1 — * 

18041.0 

17  . 

NCO+N=N2+CO 

3.30E+13 

o 

o 

o 

o 
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18. 

NCO+H=NH+CO 

5.20E+13 

o 

o 

o 

o 

19. 

NC0+0=N+C02 

8.00E+12 

0.0 

2503.0 

20. 

NC0+0=N0+C0 

4 . 2  OE+13 

o 

o 

o 

o 

21. 

NCO+M=N+CO+M 

3.30E+14 

o 

o 

54046.0 

22. 

NC0+N0=N20+C0 

6.20E+17 

-1.7 

765.0 

23. 

NC0+N0=C02+N2 

7.80E+17 

-1.7 

765.0 

24  . 

HONO+M=OH+NO+M 

2.00E+31 

-4 . 6 

51175.0 

25. 

H0N0+H=H2+N02 

2 . 00E+08 

1.6 

6616.0 

26. 

H0N0+H=H20+N0 

8.10E+06 

1.9 

3844.0 

27  . 

HONO+H=OH+HNO 

5.60E+10 

0.9 

4967.0 

28. 

H0N0+0=0H+N02 

1.70E+08 

1.5 

3030.0 

29. 

H0N0+0H=H20+N02 

1 .20E+06 

2.0 

-596.0 

30. 

H0N0+CH3=N02+CH4 

8.10E+05 

1.9 

5504.0 

31. 

H0N0+NH2=N02+NH3 

9.20E+05 

1.9 

1917.0 

32. 

HNO+M=H+NO+M 

2 . 60E+16 

0.0 

48681.0 

33. 

HN0+0H=N0+H20 

1 .30E+07 

1.9 

-953.0 

34. 

HNO+H=H2+NO 

4.50E+11 

0.7 

655.0 

35. 

HN0+H=0+NH2 

3.50E+15 

00 

o 

1 

29268.0 

36. 

HN0+0=0H+N0 

4.50E+11 

0.7 

655.0 

37. 

HNO+NH2=NH3+NO 

9.20E+05 

1.9 

-1152.0 

38. 

HN0+N02=H0N0+N0 

9.20E+05 

1.9 

-1152.0 

39. 

HN0+N0=N20+0H 

8.50E+12 

o 

o 

29586.0 

40. 

HNO+CH3=NO+CH4 

8.20E+05 

1.9 

953.0 

41. 

NO+C=CO+N 

1 . 7  OE+13 

0.0 

0.0 

42  . 

N0+C=CN+0 

1 . 1 OE+13 

0.0 

0.0 

43. 

N02+H=N0+0H 

1 .30E+14 

0.0 

357.0 

44  . 

N02+0=N0+02 

3.90E+12 

0.0 

-238.0 

45. 

N02+M=N0+0+M 

5.70E+15 

0.0 

59987.0 

46. 

N02+NH2=N20+H20 

1.50E+16 

-1.4 

268.0 

47  . 

N0+M=N+0+M 

1.40E+15 

0.0 

148428.0 

48. 

N20+M=N2+0+M 

5.70E+14 

0.0 

56093.0 

49. 

N20+0=N2+02 

1.40E+12 

o 

o 

10809.0 

50. 

N20+0=N0+N0 

2.90E+13 

o 

o 

23148.0 

51. 

NH2+N02=NH2N02 

5.60E+48 

-11.3 

11882.0 

52  . 

H0N0+H0N0=H20+N0+N02 

3.49E-01 

3.6 

12138.6 

53. 

CH30+N0=CH20+HN0 

8.40E+12 

o 

o 

2050.0 

54  . 

CH3+H02=CH30+0H 

2.28E+13 

o 

o 

o 

o 

55. 

N02+N02=N0+N0+02 

1.60E+12 

o 

o 

26123.0 

56. 

N02+N02=N03+N0 

9.60E+09 

0.7 

20900.0 

57  . 

0H+N02 (+M) =HN03 (+M) 

1.37E+15 

'D 

O 

1 

o 

o 

Low  pressure  limit: 

0 . 142  00E+2  6 

-0 .29000E+01 

0 . 00000E  +  00 

58. 

HN0+HN0=N20+H20 

4 . 00E+12 

O 

O 

5000.0 

59. 

CH3+N02=CH30+N0 

1.30E+13 

o 

o 

o 

o 

60. 

N204 ( +M) =N02+N02 (+M) 

4.05E+18 

-1.1 

12840.0 

Low  pressure  limit: 

0.19600E+29 

-0 . 38000E+01 

0 . 12840E  +  05 

61. 

N02+H0N0=N0+HN03 

6.03E+01 

o 

o 

o 

o 

62. 

N204+H20=H0N0+HN03 

2.52E+14 

o 

o 

11586.0 

63. 

H+0+M=0H+M 

6.20E+16 

T5 

O 

1 

o 

o 

64. 

H+02=0H+0 

1.59E+17 

-0.9 

16873.8 

65. 

H2+0=0H+H 

3.87E+04 

2.7 

6262.0 

66. 

H2+0H=H20+H 

2.16E+08 

1.5 

3441.7 

67. 

0H+0H=0+H20 

2.10E+08 

1.4 

-396.8 

68. 

0+0+M=02+M 

1.00E+17 

i 

i — 1 

o 

o 

o 

H2 

Enhanced  by 

2 . 900E+00 

02 

Enhanced  by 

1 . 200E+00 

N2 

Enhanced  by 

1 . 200E+00 

H20 

Enhanced  by 

1 . 850E+01 

CO 

Enhanced  by 

2 . 100E+00 

C02 

Enhanced  by 

4 . 300E+00 

H+0H+M=H20+M 

8.40E+21 

i 

N) 

O 

O 

O 

H2 

Enhanced  by 

1.900E+00 

02 

Enhanced  by 

2 . 600E+00 

N2 

Enhanced  by 

2 . 600E+00 

H20 

Enhanced  by 

9 . 500E+00 

CO 

Enhanced  by 

2 . 600E+00 

C02 

Enhanced  by 

2 . 600E+00 

H+02+M=H02+M 

7.00E+17 

00 

O 

1 

o 

o 

H02+H=0H+0H 

1.50E+14 

O 

O 

1003.8 

28 


72  . 

H02+H=H2+02 

2.50E+13 

o 

o 

693.1 

73. 

H02+H=H20+0 

5 . 00E+12 

o 

o 

1410 . 1 

74  . 

H02+0=02+0H 

2.00E+13 

o 

o 

o 

o 

75. 

H02+0H=H20+02 

6.02E+13 

o 

o 

o 

o 

76. 

H02+H02=H202+02 

4 . 2  OE+14 

o 

o 

11981.6 

77  . 

H202+M=0H+0H+M 

1 . 00E+17 

o 

o 

45411 . 1 

H2 

Enhanced  by 

2 . 900E+00 

02 

Enhanced  by 

1 . 200E+00 

N2 

Enhanced  by 

1 .200E+00 

H20 

Enhanced  by 

1 . 850E+01 

CO 

Enhanced  by 

2 . 100E+00 

C02 

Enhanced  by 

4 . 300E+00 

78. 

H202+H=H02+H2 

1.70E+12 

o 

o 

3752.4 

79. 

H202+H=H20+0H 

1.00E+13 

o 

o 

3585.1 

80. 

H202+0=H02+0H 

2.80E+13 

o 

o 

6405.3 

81. 

H202+0H=H20+H02 

7.00E+12 

o 

o 

1434.0 

82. 

CO+OH (+M) =C02+H (+M) 

2.45E-03 

3.7 

-1242.8 

Low  pressure  limit: 

0 . 117  00E  +  08 

0 . 13500E+01 

-0 . 71700E+03 

83. 

C0+H02=C02+0H 

1.50E+14 

o 

o 

23589.8 

84. 

C0+0+M=C02+M 

3.01E+14 

o 

o 

3011.5 

02 

Enhanced  by 

1 .200E+01 

N2 

Enhanced  by 

2 . 000E+00 

CO 

Enhanced  by 

3 . 000E+00 

C02 

Enhanced  by 

7 . 000E+00 

85. 

00+02=002+0 

2.50E+12 

o 

o 

47801.1 

86. 

HC0+M=C0+H+M 

1.86E+17 

o 

\ — 1 

16993.3 

H2 

Enhanced  by 

1 . 870E+00 

N2 

Enhanced  by 

9. 900E-01 

02 

Enhanced  by 

9 . 900E-01 

87. 

HOCO ( +M) =C02+H (+M) 

1.74E+12 

0.3 

32935.0 

Low  pressure  limit: 

0 . 22  900E+27 

-0 . 30200E+01 

0 . 350  62E  +  05 

00 

00 

HOCO ( +M) =C0+0H (+M) 

5 . 89E+12 

LO 

o 

33986.6 

Low  pressure  limit: 

0 . 21900E+24 

-0 . 18900E+01 

0 . 35277E+05 

89. 

CHC0+M=C0+CH+M 

6.00E+15 

o 

o 

58819.3 

90. 

CH+0=C0+H 

1 . 00E+14 

o 

o 

o 

o 

91. 

CH+0H=HC0+H 

3.00E+13 

o 

o 

o 

o 

92. 

CH+02=C0+0H 

2.00E+13 

o 

o 

o 

o 

93. 

CH2S+0=C0+H+H 

1.50E+13 

o 

o 

o 

o 

94  . 

CH2S+0=C0+H2 

1.50E+13 

o 

o 

o 

o 

95. 

CH2S+0H=CH20+H 

3.00E+13 

o 

o 

o 

o 

96. 

CH2S+02=C0+0H+H 

3.10E+13 

o 

o 

o 

o 

97  . 

CH2T+0=C0+H+H 

5 . 00E+13 

o 

o 

o 

o 

98. 

CH2T+0H=CH+H20 

1.13E+07 

2.0 

3011.5 

99. 

CH2T+0H=CH20+H 

2.50E+13 

o 

o 

o 

o 

100. 

CH2T+02=C02+H+H 

1.60E+12 

o 

o 

1003.8 

101. 

CH2T+02=CH20+0 

5.00E+13 

o 

o 

9010.5 

102  . 

CH2T+02=C02+H2 

6.90E+11 

o 

o 

501.9 

103. 

CH2T+02=C0+0H+H 

8.60E+10 

o 

o 

-501.9 

104  . 

CH2T+02=HC0+0H 

4.30E+10 

o 

o 

-501.9 

105. 

CH3+0=CH20+H 

8.43E+13 

o 

o 

o 

o 

106. 

CH3+0H=CH20+H2 

8.00E+12 

o 

o 

o 

o 

107  . 

CH3+0H=CH2T+H20 

1.13E+06 

2.1 

2437.8 

108. 

CH3+02=CH20+0H 

5.20E+13 

o 

o 

34894.8 

109. 

CH3+02=CH302 

1.70E+60 

\ — 1 

LO 
\ — 1 

1 

18785.8 

110. 

CH3+HCO=CH4+CO 

3.20E+11 

LO 

o 

o 

o 

111. 

CH30H=CH3+0H 

2 . 80E+32 

o 

LO 

1 

97992.3 

112  . 

CH4+H02=CH3+H202 

4.48E+13 

o 

o 

24641.5 

113. 

CH3+H02=CH4+02 

1.00E+12 

o 

o 

o 

o 

114  . 

CH4+0=CH3+0H 

1.90E+09 

1.4 

8675.9 

115. 

CH4+0H=CH3+H20 

1.50E+06 

2.1 

2437.8 

116. 

HC0+H=H2+C0 

7 . 2  3E+13 

o 

o 

o 

o 

117  . 

HC0+0=0H+C0 

3.00E+13 

o 

o 

o 

o 

118. 

HC0+0=H+C02 

3.00E+13 

o 

o 

o 

o 

119. 

HC0+0H=H20+C0 

1 . 0  OE+14 

o 

o 

o 

o 

120. 

HC0+02=H02+C0 

4 . 2  OE+12 

o 

o 

o 

o 

121. 

HC0+HC0=CH20+C0 

2.00E+13 

o 

o 

o 

o 

122  . 

CH20+M=HC0+H+M 

5.00E+16 

o 

o 

76481.8 

H2 

Enhanced  by 

2 . 900E+00 
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02 

Enhanced 

by 

N2 

Enhanced 

by 

H20 

Enhanced 

by 

CO 

Enhanced 

by 

C02 

Enhanced 

by 

123  .  CH20+02=HC0+H02 

124 .  CH3CHO=HCO+CH3 

125 .  HC0+CH30H=CH20+CH20H 

126.  CH20+H=HC0+H2 

127.  CH20+0=HC0+0H 

128.  CH20+0H=HC0+H20 

129.  CH20+H02=HC0+H202 

130.  CH20+CH3=HC0+CH4 

131.  CH20+M=C0+H2+M 


1 

200E+00 

1 

200E+00 

1 

850E+01 

2 

100E+00 

4 

300E+00 

1 

00E+14 

0 

0 

39914 

0 

2 

00E+15 

0 

0 

79110 

9 

9. 

.  60E+03 

2. 

.9 

13097 . 

.5 

1 . 

.  2  6E+08 

1 . 

.6 

2175  . 

.0 

3. 

.  5  OE+13 

0. 

.0 

3513. 

.4 

3. 

.  62E+05 

2  . 

.5 

-980. 

.0 

4  . 

.  00E+11 

0. 

.0 

8006. 

.7 

8. 

.  91E-13 

7  . 

.4 

-956. 

.0 

8. 

.  2  OE+15 

0. 

.0 

69550. 

.7 

H2 

02 

N2 

H20 

CO 

C02 

132.  CH30+M=CH20+H+M 

H2 

02 

N2 

H20 

CO 

C02 

133.  CH20H+M=CH20+H+M 

H2 

02 

N2 

H20 

CO 

C02 


Enhanced  by  2.900E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.850E+01 
Enhanced  by  2.100E+00 
Enhanced  by  4.300E+00 

Enhanced  by  2.900E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.850E+01 
Enhanced  by  2.100E+00 
Enhanced  by  4.300E+00 

Enhanced  by  2.900E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.200E+00 
Enhanced  by  1.850E+01 
Enhanced  by  2.100E+00 
Enhanced  by  4.300E+00 


1 . 96E+37 


1.54E+13 


0.0 


33269.6 


33054.5 


134 .  CH20H+H=CH3+0H 

135 .  CH20H+H=CH20+H2 

136 .  CH20H+0=CH20+0H 

137 .  CH20H+0H=CH20+H20 

138 .  CH20H+H02=CH20+H202 

139 .  CH20H+CH3=CH20+CH4 

140 .  CH20H+HC0=CH20+CH20 

141 .  CH20H+HC0=CH30H+C0 

142 .  CH20H+02=CH20+H02 

143 .  CH20H+H202=CH30H+H02 

144 .  CH30H+02=CH20H+H02 

145.  CH3+0H=CH30+H 

146 .  CH3+02=CH30+0 

147.  CH30+H=CH20+H2 

148.  CH30+0=CH20+0H 

149.  CH30+0H=CH20+H20 

150 .  CH 30+02 =CH20+H02 

151 .  CH30+H02=CH20+H202 

152 .  CH30+CH3=CH4+CH20 

153.  CH30+CH20=HC0+CH30H 

154 .  CH302H=CH30+0H 

155 .  CH30H+0H=CH20H+H20 

156.  CH30H+0H=CH30+H20 

157.  CH30H+H=CH20H+H2 

158.  CH30H+H=CH30+H2 

159.  CH30H+H=CH3+H20 

160.  CH30H+0=CH20H+0H 

161.  CH30H+CH30=CH20H+CH30H 

162.  CH30H+C2H3=CH20H+C2H4 

163.  CH302+H=CH30+0H 

164 .  CH302+0=CH30+02 

165.  CH302+0H=CH30H+02 

166 .  CH302+H02=CH302H+02 

167 .  CH302+HC0=CH30+H+C02 


9 

64E+13 

0 

0 

0 

0 

3 

0  OE+13 

0 

0 

0 

0 

5 

00E+12 

0 

0 

0 

0 

5 

00E+12 

0 

0 

0 

0 

1 

2  OE+13 

0 

0 

0 

0 

2 

40E+12 

0 

0 

0 

0 

1 

80E+14 

0 

0 

0 

0 

1 

20E+14 

0 

0 

0 

0 

1 

8  OE+13 

0 

0 

788 

7 

3 

01E+09 

0 

0 

2581 

3 

2 

0  OE+13 

0 

0 

44909 

2 

5 

74E+12 

-0 

2 

13934 

0 

4 

30E+13 

0 

0 

30807 

8 

2 

0  OE+13 

0 

0 

0 

0 

5 

00E+12 

0 

0 

0 

0 

5 

00E+12 

0 

0 

0 

0 

4 

28E-13 

7 

6 

-3537 

3 

1 . 

.  2  OE+13 

0. 

.0 

0 

.0 

2  . 

.  4 1E+13 

0. 

.0 

0 

.0 

1 . 

.  15E+U 

0. 

.0 

1290 

.6 

6. 

.  4  6E+14 

0. 

.0 

42997 

.1 

1. 

.  7 7E+04 

2. 

.6 

-884 

.3 

1 . 

.  7 7E+04 

2  . 

.6 

-884 

.3 

3. 

.  2  OE+13 

0. 

.0 

6094 

.6 

8. 

.  0  OE+12 

0. 

.0 

6094 

.6 

5. 

.  2  5E+12 

0. 

.0 

5329 

.8 

3. 

.  88E+05 

2  . 

.5 

3083 

.2 

1 . 

.  51E+12 

0. 

.0 

7002 

.9 

3. 

.  2  0E+01 

3. 

.2 

7170 

.2 

9. 

.  60E+13 

0. 

.0 

0 

.0 

3. 

.  60E+13 

0. 

.0 

0 

.0 

6. 

.  0  OE+13 

0. 

.0 

0 

.0 

4  . 

.  60E+10 

0. 

.0 

-2581 

.3 

3. 

.  0  OE+13 

0. 

.0 

0 

.0 

30 


168 .  CH302+CH3=CH30+CH30 

169 .  CH302+H202=CH302H+H02 

170 .  CH302+CH20=CH302H+HC0 

171 .  CH302+CH302=CH30+CH30+02 

172 .  CH 302 +CH 302 =CH30H+CH2 0+02 

173 .  CH302+CH4=CH302H+CH3 

174 .  CH302+H2=CH302H+H 

175 .  CH302+CH30H=CH302H+CH20H 

176.  CH302H+0H=CH302+H20 

177 .  CH302H+0=CH302+0H 

178.  C2H+0=C0+CH 

179.  C2H+0H=CHC0+H 

180 .  C2H+02=CHC0+0 

181 .  C2H+02=CH+C02 

182 .  C2H+02=HC0+C0 

183.  C2H2+0H=C2H+H20 

184.  CHC0+H=CH2S+C0 

185.  CHC0+0=C0+C0+H 

186 .  C2H2+0=CH2T+C0 

187.  C2H2+0=CHC0+H 

188.  C2H2+OH=HCCOH+H 

189.  C2H2+OH=CH2CO+H 

190.  C2H2+OH=CH3+CO 

191.  HCCOH+H=CH2CO+H 

192.  CH2CO+H=CHCO+H2 

193.  CH2C0+0=HC0+HC0 

194 .  CH2C0+0H=CHC0+H20 

195 .  CH2CO=CH2T+CO 

196.  C2H3+0H=C2H2+H20 

197 .  C2H3+0=CH2C0+H 

198 .  C2H3+02=CH20+HC0 

199 .  CH3CO+H=CH2CO+H2 

200 .  CH3C0+0=CH3+C02 

201 .  CH3C0+CH302=CH3+C02+CH30 

202 .  C2H4+H02=C2H3+H202 

203 .  C2H4+02=C2H3+H02 

204 .  C2H3+CH30H=C2H4+CH30 

205 .  C2H4+0=CH3+HC0 

206 .  C2H4+0=CH3C0+H 

207 .  C2H4+0=C2H3+0H 

208 .  C2H4+0H=C2H3+H20 

209 .  CH3CHO+H=CH3CO+H2 

210 .  CH3CH0+0=CH3C0+0H 

211 .  CH3CH0+0H=CH3C0+H20 

212 .  CH3CHO+CH3=CH3CO+CH4 

213 .  CH3+CH20H=C2H5+0H 

214 .  C2H5+0=CH20+CH3 

215.  C2H5+0=CH3CH0+H 

216 .  C2H5+0=C2H4+0H 

217 .  C2H5+02=C2H4+H02 

218.  C2H5+CH30H=C2H6+CH20H 

219 .  C2H6+CH302=C2H5+CH302H 

220.  C2H6+H02=C2H5+H202 

221 .  C2H6+02=C2H5+H02 

222 .  C2H6+0=C2H5+0H 

223 .  C2H6+0H=C2H5+H20 

224.  C2H5+CH30H=C2H6+CH30 

225.  NH3+M=NH2+H+M 


H20 

N2 

226.  NH3+0H=NH2+H20 

227.  NH3+0=NH2+0H 

228.  NH2+0H=NH+H20 

229.  NH2+0=NH+0H 

230.  NH2+N0=N2+H20 

231.  NH+0=N0+H 

232.  NH+0H=N+H20 

233.  NH+OH=HNO+H 


Enhanced  by 
Enhanced  by 


2.00E+13 
2.40E+12 
2.00E+12 
7 . 80E+10 
1.30E+11 
1.80E+11 
3.00E+13 
1.80E+11 
7.00E+12 
2.80E+13 
1.00E+13 
2.00E+13 
6.02E+11 
4.50E+15 
2.41E+12 
2.71E+13 
1.50E+14 
1.00E+14 
7.81E+03 
1.39E+04 
5.06E+05 
2.19E-04 
4.85E-04 
1.00E+13 
3.00E+13 
1.00E+13 
1.00E+13 
3.00E+14 
4.00E+12 
3.00E+13 
3.98E+12 
2.00E+13 
2.00E+13 
2.40E+13 
1.12E+13 
3.98E+13 
1.40E+01 
1.60E+08 
1.60E+08 
7 . 11E+08 
9.70E+06 
4 . 00E+13 
5 . 00E+12 
1.00E+13 
8.50E+10 
1.37E+14 

4 . 2  4E+13 
5.32E+13 
3.06E+13 
2.00E+12 
3.20E+01 
3.00E+11 
1.70E+13 
3.98E+13 
1.40E+00 

2 . 2  0E+07 
1.40E+01 
1.40E+16 

6 . 500E+00 

4 . 000E-01 

2.04E+06 
2.10E+13 
9.00E+07 
7.00E+12 
3.00E+20 
7.00E+13 
2 . 00E+09 
4.00E+13 


o 

o 

o 

o 

o 

o 

9942.6 

o 

o 

11663.5 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

18475 . 1 

o 

o 

26027.7 

o 

o 

13719.0 

o 

o 

1434.0 

o 

o 

6405.3 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

25095.6 

o 

o 

o 

o 

o 

o 

10492.3 

o 

o 

o 

o 

o 

o 

o 

o 

00 

CM 

501.9 

00 

CM 

501.9 

2.3 

13503.8 

4.5 

-1003.8 

4.0 

-2007.6 

o 

o 

o 

o 

o 

o 

8604.2 

o 

o 

5975.1 

o 

o 

2629.0 

o 

o 

70984.7 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

-239.0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

30425.4 

o 

o 

61496.2 

3.1 

6931.2 

1.4 

525.8 

1.4 

525.8 

1.6 

7480 . 9 

1.9 

2361.4 

o 

o 

4206.5 

o 

o 

1792.5 

o 

o 

o 

o 

o 

o 

6000.0 

o 

1 

6596.6 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

4995.2 

3.2 

9153.9 

o 

o 

14937 . 9 

o 

o 

20458.9 

o 

o 

50908.2 

4.3 

2772.5 

1.9 

1123.3 

3.1 

8938.8 

0.1 

90600.0 

2 

0 

566 

0 

0 

0 

9000 

0 

1 

5 

-457 

0 

0 

0 

0 

0 

-2  . 

.6 

924 

.0 

0. 

.0 

0 

.0 

1 . 

.2 

6 

.0 

0. 

.0 

0 

.0 

31 


234.  NH+02=N0+0H 

235.  NH+N0=N20+H 

236.  N+02=N0+0 

237.  N+OH=NO+H 

238.  N+N0=N2+0 

239.  N0+H02=N02+0H 

240.  N02=N0+0 

241.  N2H2+0=NH2+N0 

242.  N2H2+0=NNH+0H 

243.  N2H2+0H=NNH+H20 

244.  N2H2+N0=NH2+N20 

245.  NNH+0H=N2+H20 

246.  NO+NNH=N2+HNO 
2  4  7.  NH2  +NO=NNH+OH 

248.  N20+H=N2+0H 

249.  N20+0H=N2+H02 

250.  NH2+02=HN0+0H 

251.  NH+02=HN0+0 

252.  N02+0=N03 

253.  N03+0=N02+02 

254.  N03+0H=N02+H02 

255.  N03+H=N02+0H 

256 .  HN0+02=N0+H02 

257.  NH2+02=NH+H02 

258.  NH2+H02=NH3+02 

259.  NH+NO=N2+OH 

260.  NNH+0=N0+NH 

261.  NNH+0=N2+0H 

262.  NH2+N0=N20+H2 

263.  NH+0=N+0H 

264.  NH2+0=N0+H2 

265.  HCN+0=NC0+H 

266.  HCN+0=CN+0H 

267.  HCN+0=NH+C0 

268 .  CN+H20=HCN+0H 

269.  HCN+OH=HOCN+H 

270.  HCN+OH=HNCO+H 

271 .  HCN+OH=NH2+CO 

272.  CN+02=NC0+0 

273.  NC0+02=N0+C02 

274.  HNCO+CN=NCO+HCN 

275 .  CH4 (+M) =CH3+H (+M) 

Low  pressure  limit:  0.45000E+18 
TROE  centering:  0.10000E+01 

276.  CH3NH+0H=CH2NH+H20 

277 .  CH3NH+0=CH2NH+0H 

278 .  CH3NH+02=CH2NH+H02 

279.  CH3NH+0=CH30+NH 

280.  CH3NH+OH=CH4+HNO 

281 .  CH3NH+02=CH30+HN0 

282.  CH2NH+0=CH20+NH 

283.  CH2NH+0H=CH20+NH2 

284.  CH2NH+0=H2CN+0H 

285 .  H2CN+H02=CH2NH+02 

286.  CH2NH+0H=H2CN+H20 

287 .  H2CN+0=HCN+0H 

288.  H2CN+0H=HCN+H20 

289 .  H2CN+02=HCN+H02 
2  90.  H2  CN+NO=HCN+HNO 

291 .  CH3+CH3 (+M) =C2H6 (+M) 

Low  pressure  limit:  0.12600E+42 
TROE  centering:  0.62000E+00 

292.  CH2NH+M=HCN+H2+M 

293.  2H+M=H2+M 

294.  2H+H2=H2+H2 

295.  NH+N=N2+H 

296.  NH+H=N+H2 

297.  NH2+H=NH+H2 


1.00E+13 

-0.2 

4800.0 

2.40E+15 

00 

o 

1 

o 

o 

6.40E+09 

1.0 

6279.9 

3.80E+13 

o 

o 

o 

o 

3.30E+12 

0.3 

o 

o 

2 . 11E+12 

o 

o 

-479.0 

7 . 60E+18 

-1.3 

73290.0 

1.00E+13 

o 

o 

o 

o 

2.00E+13 

o 

o 

1000.0 

1.00E+13 

o 

o 

1000.0 

3.00E+12 

o 

o 

o 

o 

5.00E+13 

o 

o 

o 

o 

2.00E+12 

o 

o 

o 

o 

1.39E+12 

o 

o 

o 

o 

9.60E+13 

o 

o 

15100.0 

2.00E+12 

o 

o 

10000.0 

1.51E+12 

-3.9 

36100.0 

4 . 61E+05 

2.0 

6500.0 

1 .32E+13 

o 

o 

o 

o 

6.03E+12 

o 

o 

o 

o 

1.39E+13 

o 

o 

o 

o 

6.62E+13 

o 

o 

o 

o 

3.16E+12 

o 

o 

3000.0 

1.00E+14 

o 

o 

49997.0 

4.52E+13 

o 

o 

o 

o 

2.16E+13 

CM 

o 

1 

o 

o 

1.00E+13 

o 

o 

o 

o 

1.00E+13 

o 

o 

5000.0 

5.00E+13 

o 

o 

24481.0 

7.00E+12 

o 

o 

o 

o 

5 . 00E+12 

o 

o 

o 

o 

1.40E+04 

2.6 

4980.0 

2.70E+09 

1.6 

29200.0 

3.50E+03 

2.6 

4980.0 

8.00E+12 

o 

o 

7450.0 

5.90E+04 

2.4 

12500.0 

2.00E-03 

4.0 

1000.0 

7.80E-04 

4.0 

4000.0 

7.50E+12 

o 

o 

-389.0 

2.00E+12 

o 

o 

20000.0 

1.50E+13 

o 

o 

o 

o 

2.40E+16 

o 

o 

104913.6 

0 . 00000E+00  0 . 90806E+05 

0 . 10000E+03  0 . 13500E+04  0 . 78300E+04 


1.00E+08 

2.0 

o 

o 

1.00E+08 

2.0 

0.0 

1.00E+07 

2.0 

6300.0 

6.00E+13 

0.0 

0.0 

6.00E+12 

0.0 

0.0 

6.00E+12 

0.0 

4000.0 

1.00E+07 

2.0 

2800.0 

1.80E+05 

2.0 

14800.0 

3.16E+08 

2.0 

6100.0 

7.87E+04 

2.0 

21700.0 

1 . 00E+07 

2.0 

4000.0 

1 . 00E+07 

2.0 

6100.0 

1.00E+07 

2.0 

3700.0 

2.70E+04 

2.0 

17300.0 

1.00E+07 

2.0 

4400.0 

3.61E+13 
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7.3  Predictions  of  the  kinetics  model 


The  above  kinetics  model  is  used  here  to  predict  the  fate  of  homogeneously  premixed 
HNCO/CO2/H2O/N2O/HNO3/NO2  gaseous  compositions.  The  proposed  model  is  able  to 
simulate  ignition  at  relatively  low  temperatures.  As  shown  in  Fig.  1,  for  a  mixture  at  a 
constant  volume  and  at  an  initial  temperature  of  480  K  and  initial  pressure  of  1  atm,  an 
ignition  delay  of  19  s  is  predicted.  Here,  ignition  is  numerically  defined  as  the  time  delay  in 
the  system  at  which  point  the  temperature  profile  has  reached  a  value  within  90%  of  the 
maximum.  Typically,  detailed  chemical  kinetics  models  are  validated  against  measured 
experimental  parameters  such  as  ignition  delays,  species  profiles,  and  fundamental  flame 
speeds  that  are  obtained  under  “ideal”  laboratory  setups  with  homogeneous  mixtures  at  a 
known  temperature  and  pressure.  However,  for  hypergolic  systems,  it  is  very  difficult  to 
determine  the  precise  composition,  pressure  and  temperature  of  the  gas  phase  in  which 
ignition  is  observed  due  to  the  inherent  in-homogeneity  of  the  transiently  formed  mixtures, 
and  therefore  it  is  not  easy  to  validate  the  detailed  chemical  kinetics  model  under 
consideration.  Even  in  the  case  for  premixed  homogeneous  systems,  such  a  validation  effort 
in  many  circumstances  is  not  trivial. 

The  present  detailed  chemical  kinetics  model  to  predict  ignition  delays  in  ILs  is  first  of  its 
kind.  Although  the  predictions  are  off  by  2  orders  in  magnitude,  the  present  effort  represents 


34 


a  valuable  initial  step  towards  understanding  hypergolicity  in  ILs.  It  should  be  noted  that  the 
model  assumes  a  homogeneous  gas  mixture  in  which  ignition  occurs  isotropically  throughout 
the  system,  whereas  the  experimental  ignition  delays  reported  by  Chambreau  and  Vaghjiani  et 
al.  (2008)  are  for  transiently  formed  mixtures  that  are  highly  turbulent  as  clearly  shown  in 
their  high-speed  videos  of  the  ignition  event  which  is  highly  localized.  Refinements  to  the 
model  are  to  be  expected  in  the  forthcoming  years  as  new  experimental  data  becomes 
available  from  AFRL,  Edwards  AFB.  Similarly,  re-appraisal  of  the  rate  constant  data  and/or 
addition  of  other  pertinent  reaction  steps  as  they  become  available  in  the  literature  will 
improve  the  fidelity  of  the  model  and  enhance  its  quantitative  capabilities. 

The  detailed  chemical  kinetics  model  presented  here  can  also  be  implemented  into  a 
computational  code,  such  as  Fluent,  which  incorporates  the  physics  (fluid  dynamics)  of 
hypergolic  ignition.  However,  in  its  present  form,  our  model  probably  has  too  many  reactions 
for  its  efficient  implementation.  Therefore,  reduction  techniques  can  be  used  to  build  reduced 
chemical  kinetics  models  as  discussed  below. 
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Figure  1.  Computed  temperature  profile  against  time  for  a  mixture 
HNCO/NO2/N2O/CO2/H2O/HNO3:  1/1/2/1/1/1.  The  system  is  considered  to  be  adiabatic  and 
at  constant  volume.  P  =  1  atm,  T  =  480  K,  and  ignition  delay  is  19  s. 


7.4  Sensitivity  analysis 


Although  it  may  be  premature  to  perform  sensitivity  analysis  with  a  non-validated  detailed 
chemical  kinetics  model,  the  exercise  can  allow  for  the  building  of  a  reduced  chemical 
kinetics  model.  A  reduced  kinetics  model  is  an  abbreviated  chemical  description  of  ignition 
and  combustion.  There  are  several  approaches  for  obtaining  reduced  models  through  various 
sensitivity  analyses.  These  analyses  allow  for  the  identification  of  redundant  species  and 
reactions,  i.e.  species  and  reactions  that  do  not  significantly  affect  the  chemical  evolution  of 
the  system.  It  is  generally  recognized  that  quite  often  about  90%  of  the  initial  reactions  in  the 
detailed  model  can  be  eliminated  (depending  on  the  required  conditions  of  composition, 
temperature,  and  pressure)  to  produce  a  much  smaller  and  tractable  reduced  model  with 
fidelity  characteristics  similar  to  that  of  the  starting  model. 
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7.4.1  Brute  force  sensitivity  analyses 


These  analyses  can  be  divided  into  three  groups.  One  brute  force  sensitivity  analysis  involves 
the  running  of  the  model  by  perturbing  one  parameter  and  comparing  the  results  to  the 
predictions  of  the  unperturbed  model.  For  instance,  a  rate  constant  can  be  set  equal  to  zero  (as 
if  the  considered  reaction  was  removed  from  the  model).  Such  a  sensitivity  analysis  allows 
one  to  determine  whether  or  not  the  presence  of  a  reaction  is  important  for  ignition.  Another 
brute  force  sensitivity  analysis  uses  multipliers  for  the  pre-exponential  factor  of  a  rate  constant 
(generally  between  two  and  five).  In  this  case,  the  sensitivity  analysis  allows  one  to  assess  the 
effect  of  the  magnitude  of  the  rate  constant.  Both  techniques  are  simple  but  quickly  become 
inefficient  for  large  models  that  contain  one  hundred  or  more  reactions.  Never  the  less,  they 
are  recognized  as  the  most  utilized  approaches  for  modeling  ignition  delay  times.  A  third 
brute  force  sensitivity  analysis,  which  is  rarely  used,  requires  one  to  modify  the 
thermodynamics  data  of  key  species,  according  to  experimental  (or  computational) 
uncertainties.  This  brute  force  sensitivity  analysis  is  very  cumbersome  because  of  the  number 
of  perturbations  that  are  possible  (for  the  standard  enthalpy  of  formation  at  298.15  K,  Cp°,  and 
S°).  As  discussed  in  section  5,  we  have  begun  to  compile  the  necessary  thermochemical 
database  to  enable  us  to  perform  such  an  analysis  in  the  future. 

7.4.2  First-order  sensitivity  analyses 

These  sensitivity  analyses  have  been  developed  as  alternatives  to  brute  force  sensitivity 
analyses.  They  directly  calculate  first-order  sensitivity  coefficients  by  taking  the  derivatives 
of  the  governing  equations  and  are  particularly  suitable  for  analyzing  species  profiles. 
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However,  these  first-order  sensitivity  analyses  are  generally  less  informative  when  it  comes  to 


predicting  ignition  delays. 


7.4.3  Results 


For  the  composition  considered  in  Fig.  1,  the  present  model  suggests  that  the  presence  of  NO2 
in  the  system  has  a  significant  effect  on  the  predicted  value  of  the  ignition  delay.  With  a 
smaller  amount  of  NC>29  an  ignition  event  is  still  simulated  as  shown  in  Fig.  2  but  with  a 
longer  ignition  delay  (by  a  factor  of  about  3).  The  question  as  to  whether  or  not  NO2  is 
indispensable  for  ignition  arises.  Fig.  3  shows  that  ignition  is  possible  even  if  no  NO2  is 
initially  present.  However,  under  this  situation  the  ignition  delay  is  an  order  of  magnitude 
longer.  It  therefore  appears  that  in  this  model,  even  a  small  amount  of  NO2  promotes  ignition. 
The  present  brute  force  sensitivity  analyses  suggest  that  the  initiation  reaction,  HNO3  — >  OH 
+  N02?  may  be  an  important  process  to  consider  in  the  construction  of  a  reduced  chemical 
kinetic  model. 
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Figure  2.  Computed  temperature  profile  against  time  for  a  mixture 
HNCO/NO2/N2O/CO2/H2O/HNO3:  1/0.01/2/1/1/1.  The  system  is  adiabatic  and  at  constant 
volume.  P  =  1  atm,  T  =  480  K,  ignition  delay  is  59  s. 
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Figure  3.  Computed  temperature  profile  against  time  for  a  mixture 
HNCO/NO2/N2O/CO2/H2O/HNO3:  1/0/2/1/1/1.  No  N02  is  present  initially.  The  system  is 
adiabatic  and  at  constant  volume.  P  =  1  atm,  T  =  480  K,  ignition  delay  is  about  400  s. 
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8.  Conclusion  and  perspective 


Thermochemical  calculations  at  equilibrium  show  that  species  observed  in  the  gas  phase 
above  the  liquids  are  able  to  sustain  a  flame,  and  that  many  of  the  important  species  needed  to 
explain  hypergolic  ignition  have  probably  been  identified.  The  present  calculations  also  show 
that  the  thermal  decomposition  of  isocyanic  acid  (HNCO)  is  exothermic.  Even  if  pure  HNCO 
decomposition  cannot  sustain  a  flame,  the  exothermicity  of  the  process  most  likely  plays  an 
important  role  in  defining  the  ignition  characteristics  of  the  system. 

A  detailed  chemical  kinetics  model  has  been  proposed  to  explain  hypergolic  ignition  in  ionic 
liquid-based  systems.  The  model  is  able  to  simulate  ignition  at  relatively  low  initial 
temperatures.  Although  initial  gas  phase  presence  of  NO2  is  not  needed  to  explain  hypergolic 
ignition,  the  model  shows  that  its  presence  can  significantly  affect  the  computed  ignition 
delay  values,  and  that  introduction  of  a  very  small  amount  of  NO2  in  the  model  can 
profoundly  reduce  the  calculated  delay  times. 

A  quantitative  assessment  of  the  predictions  from  the  detailed  chemical  kinetic  model  and 
experimental  values  is  not  possible  at  this  time  since  calculations  are  done  by  considering  a 
static  adiabatic  reactor  filled  with  a  homogeneous  premixed  gas  composition  at  a  given  initial 
pressure  and  initial  temperature.  In  this  ideal  reactor,  diffusion  of  species  and  thermal  transfer 
(by  conduction  and  convection)  is  ignored.  Whereas  in  the  real  experimental  reactor,  reactive 
species  are  formed  in  the  liquid  phase  and  indeed  diffuse  from  the  hot  liquid  surface  to  the  gas 
phase.  Therefore,  gradients  of  composition,  temperature  and  pressure  exist  in  the  whole 
volume  above  the  liquid  surface,  and  the  measured  ignition  delay  will  have  contributions  from 
both  the  time  to  diffuse  (physical  delay)  and  the  time  to  react  (chemical  delay).  Hence 
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consideration  of  chemistry  alone  in  a  numerical  model  will  not  explain  quantitatively  all  of 


the  experimental  observations.  A  numerical  code  that  also  includes  diffusion  of  species,  fluid 
mechanics  and  thermal  transfer  processes  will  be  needed  to  accomplish  this.  The  detailed 
chemical  kinetic  model  presented  in  this  report  is  probably  too  big  for  its  implementation  into 
such  a  physical  code  and  therefore  the  construction  of  a  reduced  kinetics  model  will  be 
required  first.  A  complete  sensitivity  analysis  on  this  reduced  model  will  then  be  needed  to 
ascertain  the  role  of  various  oxidants  such  as  OH,  O,  NO2,  and  N2O,  etc.  in  the  system. 
Finally,  a  validation  effort  will  be  needed  to  optimize  the  reduced  kinetics  model. 
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